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Abstract 



The b — > ccs transition is usually believed to hadronize predominantly in 

— (*)— (*)— 

B -> X c Di J with the Dy> originating from the virtual W . We demonstrate 
in a variety of independent ways that other hadronization processes cannot be 
neglected. The invariant mass of cs has sizable phase-space beyond mo+rriK- 
The rate for B — ► DD KX could be significant and should not be ignored 
as was done in previous experimental analyses. We estimate the number 
of charmed hadrons per i?-decay, n c , to be ~ 1.3 to higher accuracy than 
obtained in previous investigations. Even though n c is currently measured to 
be about 1.1, observing a significant B — > DD KX would support n c ~ 1.3. 
Many testable consequences result, some of which we discuss. 
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Typeset using REVTf^X 



At present, there appears to be a conflict between experiment and theory for fitting both 
the inclusive semileptonic branching ratio and the number of charmed hadrons per B decay 

n c = 1 — B(b — > no charm) + B(b — > ccs') «1 + 5(6 — > ccs') . (1) 

The prime indicates that the corresponding Cabibbo-suppressed mode is included. Experi- 
mentally the inclusive semileptonic BR has been measured accurately to be |§ 

B(B^X£u) = (10.4 ±0.4)% , (2) 

and n c is measured as || 

n c = 1.10 ±0.06. (3) 

A value of B(b — > ccs') ~ 0.1, suggested by Eqs. ([!]) and (Q), would lead to a theoretical 
prediction of B(B — > X^i/) that is too large — i.e., inconsistent with its measured value (0). 
On the other hand, theory predicts n c ~ 1.3 when the observed semileptonic BR is used as 
input, which is demonstrated below. Thus a conflict arises between (|2|) and (|3|) 0J§]. 

Recently, Bagan et al. and Voloshin made progress on the theoretical side 0-0. Bagan 
et al. [f§-|6[] performed a complete next-to-leading order analysis of inclusive B decays, 
which included important final state mass effects in the QCD corrections. The predicted 
B(B — ► Xlv) agrees with (0), within uncertainties that are dominated by renormalization 
scale-dependences in the perturbative calculation [f§-|6[]. Simultaneously, an enhancement 
of B(b — > ccs') was found 0-0], albeit with considerable uncertainties. Table I summa- 
rizes these recent theoretical findings ||. The main sources of the large errors in those 
studies are dependence on the renormalization-scale (m&/2 < \i < 2 m^), dependence on 
the renormalization-scheme (MS versus pole mass), and uncertainties in quark masses. 
Although this theoretical analysis hints that n c may be larger than currently measured 
it is difficult to draw firm conclusions from this direct calculation of B(b — > ccs') 
in view of the large uncertainties. 
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It should also be stressed at this point, that the experimental determination of B(B — > 
Xtv) is reliable and accurate. In contrast, the measurement of n c is a sum over the inclusive 
yields of many charmed hadron species in B decays. It is thus prone to large uncertainties, 
perhaps larger than currently realized. 

Figure 1 displays the discrepancy graphically. We discuss now in some detail how the 
theoretical curve has been generated. Our objective is to draw the most accurate curve of 
n c versus semi-electronic BR with presently available theoretical calculations. We do not 
use the prediction for B{b — > ccs') because it involves large errors, but rather proceed as 
follows. We start with 

B(b — > c) = 1 — B{b — > no charm) , (4) 
where B(b — > no charm) is small, typically at the percent level. We take 

r 4 = T{b -> no charm) /r(6 -> cev) = 0.25 ± 0.10, (5) 



to account for the small fraction of b — > s + no charm JT0[ and charmless b — ► u transitions. 
Furthermore we use 

r T S = 0.25 , (6) 

T[b — > cev) 



which is in accordance with the result of Ref. [FL\ and also agrees with a recent ALEPH 
measurement [O, 



B{b -> Xrz/) = (2.75 ± 0.30 ± 0.37)% . (7) 

The last required ratio is Y(b — > cud')/T{b — > cev) where the dominant uncertainties in 
\V c b\ 2 and in fermion masses cancel. Bagan et al. ||] have presented a complete computation 
of this quantity in next-to-leading logarithmic approximation taking all final-state charm 
quark mass effects into account. Based on this perturbative calculation and also including 
nonperturbative corrections up to 0(1/ ml), the analysis of M yields, 

T(b — > cud') , . 

r ud = -prr f = 4.0 ± 0.4 . 8 

T(b —> cev) 
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Here the error comes almost entirely from the renormalization-scale uncertainty and rep- 
resents a conservative estimate when working to order 0(l/m|). Because nonperturbative 
effects at 0(1/ ml) could introduce rate-differences at the 10% level between B~ and B& 
decays governed by b — > cud ||13|| , there is considerable room for additional studies. 
Combining Eqs. (|), (|), (0), ([8]), the 6 — > ccs' branching fraction can be written as 

B(b -> ccs') = 1 - (2 + r r + r ud + ry) S(S -> X c ez/) 

= 1 - (6.50 ± 0.40) B(B -> X c e^) . (9) 

In this relation the very small contribution from 6 — ► ucs' transitions has been neglected. 
Eqs. (|^) and yield the number of charms per B decay as 

n c = 2- (2 + r T + r ud + 2ry) B(B -> X c ei/) 

= 2 - (6.75 ± 0.40) B(B -> X c ei/) , (10) 

where we note that 5(6 — *> ccs') drops out in the linear relation between n c and B — > X c ei/, 
and that the relation is largely free from uncertainties in masses of b and c quarks since the 
error is dominated by the uncertainty in r u d. Figure 1 shows the discrepancy between theory 
given by Eq. ([IT]) and experiment. 

The precisely measured semileptonic BR together with Eqs. (|9|)- ([T0|) gives 

B(b -> ccs') = 0.32 ± 0.05 , (11) 

n c = 1.30 ±0.05. (12) 

This is our central result. Our predictions for B(b — ► ccs') and for n c agree with the central 
values obtained in previous theoretical investigations |5|-[^,|2]|| but have smaller errors. As 
discussed in more detail below, such a large value of B(b — > ccs') requires a significant rate 
for B — > XX. We predict the observation of (a) B — > D^D^K modes with significant 
BWs, (b) enhanced £ + D and correlations where the primary lepton originates from one 
B and the charmed hadron from the other B in the event, and (c) enhanced DD and D D 
correlations at the T(AS) — > BB. 
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If the predicted effects will be observed, then the B(b — > ccs') is larger than currently 
determined by experiment. The measured number of charm per B will not change by those 
observations, but the larger B(b — > ccs') would indicate that the current experimental value 
of n c is underestimated. In that case, a careful re-evaluation of all errors involved in measur- 
ing n c would be in order, including re-assessments of absolute BR's of the charmed hadrons 
some of which are poorly known. On the other hand, non-observation of our predictions 
would indicate an enhancement of the b — > cud transition over the parton estimate fl4] 
and/or a larger rate than anticipated for charmless b — > s transitions fl5|]T6[ . 



Theory alone or experimental measurements alone have large inherent uncertainties for 
B{b — > ccs'). We therefore adopted a hybrid approach which uses well measured quantities 
from experiment in conjunction with reliably calculated quantities from theory to determine 
B{b — > ccs') to higher accuracy ||. 

One conventional way to determine B(b — *> ccs) is to add the inclusive yield of D s fl, [l7| . |T5 



R Ds = B(B -> D~X) + B(B -> D+X) (13) 
to the other observed final states governed by b — > ccs [IS], 



B(b -> ccs) = R Da + B(B -> E C A C X) + B(B -> (cc)X) = 

= 0.12 + 0.01 + 0.03 = 0.16 ±0.02. (14) 

(cc) denotes charmonia not seen in DDX such as J/ip, ip', r) c , r)' c , Xc, h c , 1, ^>2- Within errors, 
this agrees with the experimental measurement of n c 

B{b -> ccs') = n c - 1 + B{b -> no charm) = 0.13 ± 0.06 . (15) 

The agreement appears to support the low value of n c . 

Our determination of B(b — ► ccs') suggests a different picture as to how 6 — > ccs 
hadronizes. A systematic classification shows that five classes of hadronization can occur, 
see Table II. Conventional wisdom 0,0,0 assumes that most of the inclusive D s produc- 



tion in B decays originates from the virtual "W". Motivated by the observed inclusive 



momentum spectrum of D s in B decays and by a simple theoretical argument given below, 
we predict instead that only about 70% of the inclusive D s yield in B decays contribute to 
B — > DD~X processes. The remaining D s (about 30%) could occur in conjunction with ss 
fragmentation. We will return to this point below. 

The branching ratio for class (a) is thus depleted and becomes about 0.7 Rd s - [This 
branching ratio can be at most Rd s , which would soften our conclusion by a small amount 
only]. The branching ratios of the observed classes (a)-(c), do not add up to 30%. Thus 
class (d) must have a sizable branching fraction of about 20%, 

B(B -> DD KX) ~ 20% . (16) 

There are several interesting experimental implications. Those modes can be studied at 
CLEO and at LEP. CLEO has higher statistics, whereas LEP has the ability to separate 
one B from the other b hadron. Thus far, however, they have not been seriously searched 
for. The low Q value in this process suggests that a significant portion will be three body 



B -> D^D { * ] K . 

Because the responsible Hamiltonian is isospin zero, many isospin relations can be used to 



facilitate the observation of those modes 24 



Finally, the class (e) processes involve ss fragmentation. Their branching ratio could be 



non- negligible, at the few percent level ||22|| . A few exclusive final states would then carry 
the lion's share of the class (e) branching ratio, because of limited phase-space. 

Before proposing a number of tests, we discuss briefly a few additional indications that 
support our hypothesis from 

(a) a naive Dalitz plot analysis p5| , 



(b) measured inclusive kaon yields in B decays, and 

(c) measured inclusive D momentum spectra in B decays. 

Figure 2 shows the b — > ccs Dalitz plot resulting from the (V — A) x (V — A) matrix 
element, where the initial and final spins were averaged and summed. In this simple model, 



the cs system hadronizes as a D s X for invariant cs masses below mjj + vhk- In contrast, 
for 

rrics > m D + m K 

the cs is not seen as a D~X but rather as D KX. The Dalitz plot region contributing to 
D s production in b — > ccs decay is m 5s < mo + tuk, and one obtains 

r( ?r + -f - 0-35 . 

1 (o — > ccs) 

This argument suggests that a large fraction of b — > ccs transitions has not been ac- 
counted for in previous investigations P/XTf. (See however the analyses of Refs. p5|j26 



which reach similar conclusions to ours.) Of course, the naive Dalitz plot argument is rather 
crude. It does not address issues of hadronization, resonance bands and their interferences, 
QCD-corrections, and interferences between penguin-amplitudes (b — ► s) with the dominant 
spectator-amplitude (b — > ccs). Nevertheless, the Dalitz plot conveys the important message 
that a significant fraction of b —>■ ccs processes could be seen in DD KX. 

The surplus of the inclusive kaon yield in B decays beyond all the conventional sources 
again indicates a significant B(B — ► DD KX) fl22f . The indication is further strengthened 



by the large observed i^-flavor correlation with its parent I?-flavor at time of decay p7| , |28 
The flavor of the kaon in B — » DD KX is 100% correlated with its parent 6-flavor. The mo- 
mentum spectra of the inclusive D yields in B decays indicates an excess of low momentum 
D J s over conventional sources ||29|| . A natural explanation can be found in B — > DD KX. 



We are now ready to suggest several tests. In addition to the "indirect" measurement 
using B(b —> ccs') ~ n c — 1 which involves large errors, we suggest to directly determine 
B{b — > ccs') by adding up the "wrong-sign" charm in tagged B decays H§, 

B(b -> ccs') w B(b -> c) = B(B -> D S ~X) + B(B -> M) + B(B -> A C X) + 

+ B(B -> SIX) + B(B -> (cc)X) . (17) 

The traditional lepton and K ± tags could be supplemented by other tags, such as K* and jet 
charge techniques. Further, the number of DD and DD S events per T(4S f ) — > 5-B decay can 



be combined with the single, inclusive D and D s yields in untagged B decay to determine 
B(B — > DX) and B(B — > D~X) [|J. Of course, mixing; effects must be corrected 

for PBfl . No tagging is required to measure B(B — > (cc)X). 

A sizable B(B — > DD KX) would show up as a D^K (from cs) enhancement. The 
background at the T(45) is much reduced because 

T(45) — > ~BB —>T)^K 

(18) 

which naturally yields a .D-fT correlation, while its DK correlation is suppressed. The Dalitz 
plot allows to enhance the D^K signal correlation further by assuming 

dT dT 

dm D(*)K dm2 cs ' 

The invariant mass spectrum of the b — > ccs transition indicates that D^K (from cs) tends 
to have a large invariant mass, see Fig. 2. 

The inclusive D s yield in B decays, Rd s , has two roughly equal contributions. Figure 



3 shows the measured momentum spectrum |Tj|. Whereas the high peak is dominated by 



the exclusive two-body modes B — > £)(*)£)(*) ; the underlying dynamics of the remainder 
had been unclear. The factorization assumption is successful in predicting ratios of rates for 
the two-body modes B -> DWdW" 0. Thus we assume factorization and predict tha'; 
b — » c + -D^*- 1- is dominated by the exclusive two-body decays 5 — > £)(*) £)(*)" in analogy to 
semileptonic decay of i? mesons. We calculate that 

^V^ = 0.7 ±0.2, 19 

r(^c + #) 

where the quoted error refers to a variation in the 6-quark mass, 4.4 < mj < 5.2 GeV, and 
in the slope of the Isgur-Wise function |[30| , p 2 = 0.84 ± 0.15. The numerator is the sum 
over the four exclusive two-body rates obtained [pT] , p2| by using the heavy quark limit [|33 



The denominator is the sum of two rates b — > c + D s and b — > c + D*. It treats the b — ► c 
transition as if it were that of free quarks |15| . The decay constant fo a , the CKM elements 
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and the factorization parameter j34[] a± all cancel in the ratio. The prediction Eq. (|T9| ) can 



currently be tested since the ratio T(B £>M£>W-)/r(& -> c + D^~) is an observable |35|] 



in which the uncertainty due to B(D S — > 07r) cancels. The prediction Eq. (|19|) together with 



the measured ratio 19 



B(B -> DWDM") 

= 0.46 ± 0.04 , (20) 



i? 



yields that [0 



5(5 -> 55;X) « 5(6 -> c + DM") = (0.7 ± 0.2)R Ds . (21) 

The remainder of the inclusive D s yield in 5 decays [Rd s — 5(5 — > DD~X) = (0.3±0.2)5dJ 
could be a significant fraction of the lower momentum D s mesons. One sizable source for it 



could be the 6 — *> ccs transition with ss fragmentation |22 



B{b -> ccs + ss) « 0.01 - 0.03 . (22) 

One generally expects one D~ per such a transition, as long as D**~ and higher D~ resonance 
production in b — > ccs + ss transitions is negligible. The total D~ production in flavor-tagged 
B decays is thus expected to be 

5(5 -> D~X) « 5(6 -> ccs + ss) + 5(5 -> 55 S 7X) « 0.1 . (23) 

The 5+ yield in flavor-tagged 5 decays is governed by the 6 — > c transition with ss frag- 
mentation, and may be non-negligible 

5(5 -> 5+X) = 5 Ds - 5(5 -> 5;X) ~ 10~ 2 . (24) 

For a model of the relative contributions to the 5+ yield from 6 — > cud, civ, ccs transitions 
with ss fragmentation, we refer the reader to Ref. |[22|| . The 5+ yield in flavor-tagged 5 
decays has been traditionally neglected [0,[Lj],[J. 

In conclusion, by combining reliable theoretical calculations and precise experimental 
measurements ||, we obtain a more accurate estimate of 5(6 — > ccs') and of n c than 
previous investigations P,p7|j2]j3|j5H7| ■ We predict 
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B(b -> ccs 1 ) = 0.32 ± 0.05 and n c = 1.30 ± 0.05, (25) 

which is significantly larger than the low experimental value n c \ exp = 1.10±0.06. We believe 
that (|25| ) is on firm ground, and expect an increase in the measured n c in the future. Our 
prediction can be tested in a variety of ways. First we advocate to measure B(b — > ccs') by 
counting up the number of anticharmed hadrons (the "wrong" charm flavor) per 5-decay. 
A sizable BR for B — > DD KX is our second prediction. It shows up as a large l~ D 
and i + D correlation after removing B°-W mixing effects [pql , where the primary lepton 



comes from one B hadron and the charmed meson from the other i?-hadron in the event. 
It can also be seen by observing the exclusive modes B — > D^D* K, and/or by searching 
for D^K (from cs) correlations. There are many additional implications, consequences and 
tests which we hope to discuss in a forthcoming report |22| . 
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TABLES 

TABLE I. The predicted semileptonic branching ratio, the B(b — > ccs') and n c taken from 
Bagan et al. [6]. 



Scheme 


B(B -> X c lv) 


B(b -> ccs') 


n c 


MS 


0.112 ± 0.017 


0.35 ± 0.19 


1.35 ± 0.19 


Pole mass 


0.120 ± 0.014 


0.27 ± 0.07 


1.27 ± 0.07 



TABLE II. The five classes of hadronization of b — > ccs. (cc) denotes charmonia not seen in 
DDX, and class (e) involves ss fragmentation. 



Class 


Mode 


Si? 


Reference 


(a) 
(b) 
(c) 
(d) 
(e) 


B DD-X 
B ~ C A C X 
B -> (cc)KX 
B -> DD KX 
b — > ccs + ss 


(0.7±0.2)i2 Ds w 0.08 
0.01 
0.03 
~ 0.2 
~ /ew x 10 -2 


[21| 

in 


Total: 


6 — > ccs 


0.31 ± 0.05 
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FIGURES 
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FIG. 1. Number of charm per B decay (n c ) is plotted against the B meson semileptonic branch- 
ing ratio. The uncertainty in the theoretical prediction is indicated by dashed lines. 
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FIG. 2. Dalitz plot of the decay b — > ccs as a function of u{= rn 2 cs /mfy and s(= m| s /m^). The 
projection onto the s axis is shown at the bottom where the D K threshold is indicated by an 
arrow. 
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FIG. 3. Momentum spectrum of inclusive D s mesons produced in untagged B decays at the 
T(45) as measured by the CLEO collaboration. The parameter x is defined by x = PD S /Pmax 
where p^ax = -^bcam ~ ■ The continuum background has been subtracted. 
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